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Electronic structure of graphene functionalized with boron and nitrogen 
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We present a theoretical study of the structural and electronic properties of graphene monolayer functionalized 
with boron and nitrogen atoms substituting carbon atoms. Our study is based on the ab initio calculations 
in the framework of the density functional theory. We calculate the binding energies of the functionalized 
systems, changes in the morphology caused by functionalization, and further the band gap energy as a 
function of the concentration of dopants. Moreover, we address the problem of possible clustering of dopants 
at a given concentration. We define the clustering parameter to quantify the dependence of the properties 
of the functionalized systems on the distribution of B/N atoms. We show that clustering of B/N atoms in 
graphene is energetically unfavorable in comparison to the homogenous distribution of dopants. For most of 
the structures, we observe a nonzero energy gap that is only slightly dependent on the concentration of the 
substituent atoms. 

Keywords: functionalized graphene, boron and nitrogen, density functional theory, energy band gap 



I. INTRODUCTION 



Graphene, a two-dimensional monolayer formed out of 
sp^ hybridizated carbon atoms ordered in the hexago- 
nal lattice, is a material with unique properties. Due 
to the hexagonal symmetry, its valence and conduction 
bands with linear dispersion cross at the so called K point 
and determine semimetallic character of graphene and its 
extremely high electron mobility.^ This, in connection 
with excellent mechanical properties, renders graphene 
an ideal candidate for applications in fiexible electronics. 
However, the ability of generating controllable band gap 
is a prerequisite for effective applications of graphene in 
transistor based electronic devices. Therefore, an effec- 
tive functionalization that would open the zero energy 
band gap in pristine graphene without significant de- 
terioration of the remaining advantageous properties is 
searched for, and practically any thinkable way of reach- 
ing this goal is currently investigated. 

There exist three general methods of opening the band 
gap of graphene: application of external electric field, us- 
age of graphene nanoribbons, and functionalization. One 
of the practical realizations of the third technique consid- 
ered in this work is doping the graphene layer with boron 
and nitrogen^^. These elements, standing nearby carbon 
in the periodic table, should act as electron acceptors (B) 
and donors (N), allowing for fabrication of contemporary 
devices in the manner of the CMOS silicon technology^. 
Here, we report new results of our theoretical studies of 
electronic properties of graphene doped with boron and 
nitrogen atoms. 
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II. CALCULATION DETAILS 

Our studies of functionalized graphene layer are 
based on ah initio calculations within the framework of 
the spin polarized denstity functional theory (DFT}^i^. 
Generalized Gradient Approximation of the exchange- 
correlation functional in Perdew-Burke-Ernzerhof (PBE) 
parametrization has been applied^^. Calculations with 
periodic boundary conditions have been performed us- 
ing Siesta packageiiii^. Valence electrons have been rep- 
resented with double zeta basis sets of orbitals local- 
ized on atoms; with polarization functions also included. 
The infiuence of core electrons has been covered within 
pseudopotential formalism. Norm-conserving Troullier- 
Martins nonlocal pseudopotential^ used in this study 
have been cast into the Kleinman-Bylander^^ separable 
form. The energy cut-off determining the density of the 
utilized real space grid has been set to 800 Ry.The Bril- 
louin zone has been sampled in the 5x5x1 Mokhorst-Park 
scheme. For the band structure and density of states cal- 
culations k-sampling changed to 15x15x1 scheme. Cal- 
culations were performed within the supercell geome- 
try with graphene layers separated by a distance large 
enough to eliminate any interactions. Structural opti- 
mization has been conducted using the conjugate gra- 
dient algorithm to achieve residual forces acting on the 
atoms lower than 0.001 eV/A. 

We have performed calculations for 5x5 supercells con- 
taining 25 graphene's primitive unit cells (i.e., consisting 
of 50 atoms). Such supercell has been chosen in order to 
examine wide range of dopant concentrations and a vari- 
ety of possible distributions of the substituent atoms, and 
also to avoid band gap reduction caused by certain super- 
cell symmetries'^. In the described graphene supercell, 
one to ten B/N atoms have been introduced (leading to 
the corresponding concentration of 2-20%). In the case of 
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two substituent atoms in a supercell all the 11 symmetri- 
cally nonequivalent configurations of the atoms have been 
examined. In the remaining cases 12 different randomly 
chosen configurations have been taken into account. Op- 
timized geometries and band structures, as well as re- 
lated properties - energy band gaps, binding energies, 
and shifts of the Fermi level with respect to the top of 
the valence band and bottom of the conduction band for 
B and N dopants, respectively, have been also obtained 
for all the considered concentrations and symmetrically 
nonequivalent configurations of substituent atoms. All 
the values, for each concentrations, have been averaged 
over the investigated configurations and compared with 
relevant results for B/N substitution. 

The binding energy per atom, which has been calcu- 
lated according to the formula: 
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is used as the measure of the stability of the systems 
studied. Etot is total energy of functionalized graphene 
and Eatom,a is the total energy of the free atom of type 
a (i.e., C, N, or B). 

To quantify the effect of B/N atoms distribution over 
the graphene lattice on the calculated properties of the 
doped systems, a special parameter measuring the level of 
clustering (further called clustering parameter) of atoms 
has been introduced. This parameter is, in general, lin- 
early proportional to the sum of n-1 shortest nonequiva- 
lent distances between B/N atoms (where n is the number 
of B/N atoms in a supercell) and it is normalized to reach 
the value of -1 for maximal and 1 for minimal clustering 
in a system with geometry of pristine graphene. The 
described clustering parameter is given by the following 
mathematical formula: 



c = 



2y-y 



En 
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where J]] is a sum of n-1 shortest nonequivalent dis- 
tances between B or N atoms in the system measured 
after geometry optimization, Yl,max Smin max- 
imal and minimal value that ^ reaches for a given n 
for ideal hexagonal lattice. We have assumed that both 
substituent atoms, B and N, create bonds among them- 
selves one another equal or greater than C atoms in ideal 
hexagonal lattice and change the symmetry. Defined in 
this way, the clustering parameter for two dopants in 
the supercell (n=2) is the most negative, if B/N atoms 
are linked by a chemical bond, and the most positive if 
they are the furthest apart as possible for a given su- 
percell. For example for the case of two N atoms, in 
considered 5x5 supercell, in which the the shortest dis- 
tance in minimal clustering is 7.143 A (J^) after geometry 
optimization, the minimal and maximal distance in pris- 
tine graphene are equal to 1.433 A (Xlmin) 7.165 A 
(Z^max)' respectively; the clustering parameter c is 0.99. 



III. RESULTS AND DISCUSSION 

We start the presentation of results by considering 
structural properties of functionalized graphene mono- 
layers. The geometries of optimized structures for two 
boron (c=0.97) and nitrogen (c=0.99) atoms in 5x5 
graphene supercell, the cortesponding electron density, as 
well as bond lengths and angles between the substituent 
atom and carbon atoms are presented in Fig. [TJa) and 

(b). 

Boron atoms in comparison to nitrogen atoms have a 
greater influence on the geometry of the doped graphene 
layer, mostly due to the fact that the covalent radius of 
boron is larger than that of a carbon atom, whereas nitro- 
gen has the covalent radius similar to carbon. It turns out 
that the graphene monolayer functionalized with boron 
atoms is no longer flat because B atoms stick out from the 
surface. Similar observation was pointed out in Refs^^i^*^. 

Analyzing isolines in Fig.[T]^a) and (b), one can see re- 
gions of electron density depletion around the positions 
of boron atoms, as well as an increase of electron den- 
sity in the vicinity of nitrogen atoms. Substituent atoms 
influence only the electron density around those carbon 
atoms with which they are linked by a chemical bond. 
Electronic density around the second and further neigh- 
bors remains unchanged. 

To measure stability of the considered structures, we 
have plotted, in Fig. [U^c), the dependence of the bind- 
ing energy per atom on concentration of B/N atoms. 
Both types of structures do not segregate after the geom- 
etry optimization and binding energies are more negative, 
for B-functionalized graphene than for N-functionalized 
one. Therefore the boron doped structures are more 
stable and the differences in binding energies between 
B- and N-doped graphene increase with growing con- 
centration of dopants, from 0.02 eV up to 0.25 eV for 
2% and 20% concentration, respectively. Our observa- 
tions are in agreement with previous theoretical work^^. 
which shows that substitution with boron atoms costs 
less energy than substitution with nitrogen atoms. Mo- 
rover, boron doped monolayer graphene and nitrogen 
doped monolayer graphene were already sythesized by 
CVD method^iiiii and are stable enough for studying 
their physical and chemical properties. 

To account for the influence of the distribution of B/N 
atoms in a supercell on stability of the functionalized 
graphene, the clustering parameter c has been calculated 
according to the formula given by Eq. [2] for each concen- 
tration of dopants. Exemplary results for 12% concentra- 
tion both for boron and nitrogen, are shown in Fig. [2l As 
one can see, from the (dotted) lines fitted to each case, 
the binding energy barely dependends on the clustering 
level. For nitrogen, as well as for boron, the configuration 
of dopants being closer to its own species is less energeti- 
cally preferable than being surrounded and separated by 
carbon atoms. To have a better insight into this prob- 
lem, further calculations including more symmetrically 
nonequivalent realizations of each concentration should 
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Number of doped atoms per (5x5) supercell 

FIG. 1. The optimized structure and electron density de- 
picted for clustering parameter close to 1 for two B/N atoms 
in 5x5 graphene supercell. Depletion of electron density in 
the vicinity of B atoms (a) and its higher values around the 
position of N atoms (b) can be observed. Bond lengths and 
angles between B/N and C atoms are also denoted for each 
case, (c) The dependence of the maximum, minimum and av- 
eraged binding energy per atom (with standard deviations) for 
functionalized graphene on the concentration of B/N atoms. 



functionalized graphene 
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FIG. 2. The dependence of binding energy per atom on clus- 
tering parameter for 12% concentration of B and N atoms. 



ing of the band gap is observed in case of the majority 
of configurations, with maximal gap amounting to al- 
most 0.5 eV for B- and 0.6 eV for N-doping. The mean 
value of the band gap slightly increases with the growing 
number of substituent atoms. This trend is more pro- 
nounced for lower concentrations, where the band gaps 
range from about 0.1 eV for two atoms per supercell to 
nearly 0.25 eV for 7-10 atoms. The average band gap 
saturates slightly for higher concentrations. In the case 
of six B/N atoms in the supercell, a decrease in an av- 
erage energy gap width by over 0.5 eV can be noticed 
in comparison with five B/N atoms. It can indicate that 
averaging over 12 random configurations may not be suf- 
ficient to eliminate the strong influence of some particu- 
lar configurations on the average values of the band gap. 
The results of our calculations signalize a limited possi- 
bility of band gap engineering by means of increasing the 
concentration of functionalizing species. 



IV. CONCLUSIONS 



be supportive. 

Fermi level shift with respect to VBT*^ /CBB*^ result- 
ing from B/N functionalization and loss/gain of one elec- 
tron per substituent atom is depicted in Fig. [3l Similar 
observation, but for smaller concentrations was denoted 
in Refs.'^'^^i^ 

Within each of the considered concentrations of B/N 
atoms, we have performed electronic structure calcula- 
tions for different configurations of dopants correspond- 
ing with this concentration. In the next step, we have 
performed averaging of the calculated band gaps. The 
average band gaps for N-doped graphene together with 
the maximal and minimal values of the band gaps are 
presented in Fig. |4j For B-doped graphene, the picture 
is analogous and is therefore not presented here. Open- 



In this paper we report the results of ab initio calcu- 
lations for graphene functionalized with boron and ni- 
trogen atoms. We investigate structures with one to 
ten borons or nitrogens in a graphene supercell for vari- 
ous symmetrically nonequivalent configurations. All the 
structures are stable with slightly more negative binding 
energies in the case of B-doping. Doping with B atoms 
introduces more changes in the morphology of the func- 
tionalized structures than doping with N atoms. Both 
types of dopants prefer to be homogenously distributed 
over the lattice rather than clustered. We conclude that 
the effects of the band gap oppening in the boron and 
nitrogen functionalized graphene are qualitatively and 
quantitatively very similar. However, the magnitude of 
the band gap could be strongly dependent on a particular 
configuration of dopants. To clarify this issue, the role of 
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FIG. 3. The maximum, minimum and averaged (with stan- 
dard deviations) shift of the Fermi level from (a) the top of 
the valence band (VBT) for B and (b) the bottom of the 
conduction band (CBB) for N in function of concentration 
of substituent atoms. Inset: band structure and density of 
states for the case of one substituent atom in 5x5 supercell. 
Shifting of Fermi energy is marked with red arrows. 



disorder should be further studied. 
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